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PREDICTING  THE  PERFORMANCE 
OF  SOLAR  ENERGY  SYSTEMS 


1 INTRODUCTION 


Background 

The  technical  feasibility  of  heating  and  cooling 
buildings  using  flat-plate  solar  collectors  has  been 
established  both  in  theory  and  practice.  Although 
improvements  in  component  and  system  design  are 
forthcoming,  one  can  approach  solar  heating  and 
cooling  technology  with  full  confidence  that  a practi- 
cal. reliable  system  can  be  constructed.  The  con- 
struction phase  of  solar  systems  requires  little  more 
skill  than  is  presently  required  to  install  conventional 
heating  and  cooling  systems;  the  design  phase,  how- 
ever. is  considerably  more  complex. 

Sunlight  is  inherently  an  intermittent  source  of 
energy  at  the  earth’s  surface.  For  this  reason,  solar 
energy  systems  for  heating  domestic  hot  water  or  for 
heating  and  cooling  buildings  are  not  and  frequently 
cannot  be  designed  to  meet  the  full  demands  of  the 
building  or  buildings  being  served.  System  design  is 
therefore  unconventional  in  that  solar  energy  sup- 
plies some  but  not  all  of  the  building’s  energy  re- 
quirements. The  fundamental  problem  in  analyzing 
solar  energy  systems  is  determination  of  the  collector 
array  area  which  will  provide  the  greatest  cost  bene- 
fit in  meeting  the  highest  fraction  of  the  building's 
energy  load.  In  the  past,  because  of  the  uncontrolled 
nature  of  solar  energy,  extensive  and  costly  computer 
studies  (taking  into  account  hourly  weather  data 
from  the  site  in  question)  have  been  required  to 
evaluate  proposed  solar  projects  and  to  design  solar 
energy  systems.  The  tools  used  for  these  studies  have 
generally  been  proprietary  and  thus  not  available  to 
Corps  of  Engineers  Districts.  Those  which  have  been 
available  are  extremely  difficult  to  use  and  are  not 
ready  for  general  use  by  design  engineers. 

Hence,  there  is  a basic  need  for  a simple,  con- 
venient method  for  making  the  necessary  design 
calculations  and  for  assessing  energy  and  life-cycle 
costs.  In  particular,  simplified  techniques  for 
making  a preliminary  energy  and  life-cycle  cost 
analysis  are  needed  to  minimize  the  number  of 
expensive  computer  studies  required.  Where  com- 
puter studies  are  required,  user-oriented  computer 
programs  are  necessary. 
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Objective 

this  work  is  part  of  a study  to  develop  a straight- 
forward procedure  to  allow  District  Engineers  to 
perform  an  energy  and  life-cycle  cost  analysis  and 
final  design  calculations  for  individual  solar  energy 
system  applications.  The  procedure  was  to  permit 
(I)  preliminary  solar  energy  system  performance 
assessment  Irom  charts  and  graphs  and  (2)  the  com- 
puter-aided evaluations  of  solar  energy  system  per- 
formance necessary  for  final  design,  and  was  to  be 
applicable  to  water  and  air  heating  systems  and 
water  heating  and  cooling  systems.  The  purpose  of 
this  report  is  to  describe  how  to  use  the  graphical 
and  computer  tools  developed  in  response  to  the 
study  objective. 

Approach 

The  approach  used  in  this  research  consisted  of 
the  following  steps: 

1.  Develop  a computer  simulation  program  for 
analyzing  the  performance  of  solar  energy  systems. 

2.  Use  the  solar  simulation  program  to  determine 
the  performance  of  typical  solar  systems  In  meeting 
predetermined  hourly  building  loads  at  each  of  five 
sites  in  the  continental  United  States  using  hourly- 
incident  radiation  data. 

3.  Study  the  effects  on  solar  system  performance 
of  varying  solar  system  design  parameters  (including 
collector  area,  collector  tilt  angle,  storage  tank 
volume,  and  heat  exchanger  effectiveness). 

4.  Based  on  the  results  of  the  parametric  studies, 
develop  a graphical  method  for  estimating  solar 
system  performance.1 

5.  Develop  a user-orien'ed  computer  simulation 
program  for  making  rigorous  solar  system  perform- 
ance calculations  during  the  final  design  phases.1 

b.  Develop  recommended  procedures  for  using 
the  graphical  and  computer  simulation  techniques 

‘D.  C.  Hittle.  D.  F.  Holshouser,  and  G.  N.  Walton.  Interim 
Feasibility  Assessment  Guidance  for  Solar  Heating  and  Cwding  of 
Army  Buildings.  Technical  Report  E-91  ADA026588  (U.S.  Army 
Construction  Engineering  Research  Laboratory  (CERL).  1976). 

2 Building  Loads  Analysis  and  System  7 hermttdynamics 
(BLAST)  Program  Users  Manual.  Draft  Technical  Report 
(CERL.  1977);  and  Building  Loads  Analysis  and  System  Thermo- 
dynamics ( BLAST)  Program  Reference  Manual.  Draft  Technical 
Report  (CERL.  1977). 


for  making  solar  system  energy  and  life-cycle  cost 
analysis,  and  for  optimizing  solar  system  designs 
based  on  lowest  system  life-cycle  cost.  This  report 
describes  these  procedures. 

Organization  of  Report 


2.  Annual  and  monthly  energy  requirements,  Ql. 
For  the  purposes  of  the  universal  curve,  the  energy 
requirements  are  defined  as  the  total  thermal  energy 
in  the  form  of  hot  water  required  by  the  building's 
domestic  hot  water  heating  system,  space  heating 
system,  or  absorption  chiller. 


Chapter  2 introduces  the  graphical  method  for 
computing  expected  performance  of  solar  collector 
arrays,  and  Chapter  .1  presents  examples  of  its  use. 
Chapter  4 describes  a method  for  optimizing  the 
solar  energy  system  on  the  basis  of  life-cycle  cost. 
Chapter  5 describes  the  computer-aided  method  for 
solar  energy  systems  w hich  is  recommended  for  per- 
forming design  calculations  on  larger  projects. 
Chapter  6 outlines  some  practical  considerations 
relevant  to  the  application  of  solar  technology. 

THE  UNIVERSAL  SOLAR  SYSTEM 
2 PERFORMANCE  CURVE— 

A GRAPHICAL  APPROACH 


3.  Percentage  of  energy  requirements  met  by  solar 
energy,  p.  T his  term  is  the  percentage  of  Qj_  met  by 
the  solar  energy  system. 

Figure  1 shows  the  form  of  the  three  universal 
curves  presented  in  this  report — one  for  domestic  hot 
water  heating  only  (shown  in  Figure  1).  one  for  heat- 
ing only  (with  or  without  domestic  hot  water  heat- 
ing). and  one  for  heating  and  cooling. 

The  X axis  is  the  solar  system  performance  param- 
eter. Ps,  which  is  defined  as  the  ratio  of  annual  or 
monthly  incident  radiation  on  the  collector  array  to 
the  annual  or  monthly  energy  requirements  of  the 
building,  as  given  by  Eq  1: 


Introduction 

During  the  development  of  a solar  system  simula- 
tion model  and  design  methodology,  the  U.S.  Army 
Construction  Engineering  Research  Laboratory 
(CERL)  performed  several  hundred  solar  system 
simulations  for  typical  Army  buildings  in  various 
parts  of  the  country.  Analysis  of  the  solar  system 
performance  curves  for  these  systems  indicated  that 
with  proper  normalization,  the  performance  of  a 
given  solar  system  for  the  various  buildings  in  all 
locations  could  be  represented  by  a single  universal 
performance  curve. 

Use  of  this  curve  requires  calculation  of  building 
energy  loads  and  compilation  of  on-site  weather 
data.  Once  this  information  is  available,  the  curve 
relates  solar  system  performance  to  the  collector 
array  area.  In  order  to  apply  the  curve,  however,  an 
understanding  of  the  following  terms  is  required: 

1.  Annual  or  monthly  incident  solar  radiation 

flux  density.  Qc.  This  is  the  solar  flux  density  on  the 

tilted  collector  array  in  Btu/square  foot/month  or 

Btu/square  foot/year  (langleys/month  or  langleys/ 
year).  Note  that  solar  data  for  a particular  site  are 

given  in  terms  of  horizontal  radiation  densities. 

Thus,  these  numbers  must  be  corrected  for  the  tilt 

angle  of  the  collector  array.  The  method  for  making 

these  corrections  is  presented  in  the  Solar  Radiation 

Data  section  of  this  chapter. 


P 


s 


0CAC 

"or 


where  Ac  = collector  area. 


I Eq  1 ) 


The  Y ordinate,  q,  is  the  fraction  of  Ql  supplied  by 
the  collector  system.  For  example,  a 200-man  bar- 
racks in  Topeka.  KS  requires  1.73  x 101 * * * * * * * 9  Btu  (1.83  x 
109  kJ)  of  energy  annually  to  heat  hot  water.  Use  of 
the  universal  curve  permits  determination  of  the 
percentage  of  the  load  which  can  be  met  by  2400  sq 
ft  (223  m2)  of  collectors  tilted  at  39  degrees. 


The  solar  radiation  flux  density.  Qc,  on  a collector 
tilted  at  39  degrees  is  6.0  x 10s  Btu/sq  ft  (5.9  x 104 
kJ/m2).  (The  method  for  obtaining  solar  flux  density 
on  a tilted  surface  is  described  in  the  Solar  Radiation 
Data  section.)  From  Eq  1,  Ps  = .83;  from  Figure  1. 
q is  .48.  Thus,  2400  sq  ft  (223  m2)  of  collector  would 
meet  48  percent  of  the  load;  i.e.,  the  collector  array 
would  supply  (.48)  (1.73  x (09  Btu)  = 8.3  x 1 0®  Btu 
(8.8  x 10*  kJ),  and  the  auxiliary  system  would  have  to 
supply  9 x 10*  Btu  (9.5  x 1 0®  kJ)  with  a conventional 
hot  water  heating  system. 


If  the  efficiency  of  the  auxiliary  energy  supply 
system  (an  oil  boiler,  for  example)  is  80  percent,  the 
solar  energy  system  would  save  10.4  x 109  Btu 
(10.97  x 109  kJ)  of  fuel  energy.  If  the  heating  value  of 
oil  is  150.000  Btu /gal  (41  809  kJ/f)  at  a cost  of 
$0. 40/gal  ($0. 1 1 /(),  the  annual  dollar  savings  would 
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I'igure  1.  Universal  curve  for  hot  water  heating  with  solar  energy. 


he  $27. ’(X).  This  savings  could  be  used  to  offset  the 
capital  cost  of  the  solar  system.  With  this  informa- 
tion, a life-cycle  cost  analysis  for  this  candidate  solar 
system  can  be  made.  The  procedure  can  be  repeated 
for  different  collector  areas  until  the  collector  area 
giving  the  maximum  cost  benefit  is  found. 

It  should  be  pointed  out  that  the  universal  curve 
contains  several  assumptions.  First,  a single-cover 
flat-plate  collector  with  an  emissivitv  of  0.10  and  an 
absorptivity  of  0.90  was  used  as  the  reference  collec- 
tor to  develop  the  universal  curve.  (A  method  for 
adjusting  the  collector  area  for  collector  types  other 
than  the  reference  collector  is  presented  in  Chap- 
ter 3.)  Second,  the  storage  tank  was  assumed  to  hold 
16  lb  of  water  per  square  foot  of  collector  (78  kg/m2 
or  about  2 gal/sq  ft).  This  was  found  to  be  near  opti- 
mum in  CERL  simulation  studies.  Variations  in 
tank  volume  can  be  evaluated  in  the  final  design 
phase  using  the  computer  simulation  model  (see 
Chapter  5). 

Obtaining  Required  Input  Information 


energy  demand.  Ql>  and  the  radiation  flux  density. 
Qc.  The  following  sections  describe  how  this  infor- 
mation can  be  obtained. 

Energy  Consumption  Data 

The  type  of  data  required  to  determine  and  the 
sources  of  that  data  depend  on  the  type  of  system 
being  considered  Since  solar  energy  systems  pro- 
duce thermal  energy  in  the  form  of  hot  water  (or  hot 
air  in  the  case  of  hot  air  solar  heating  systems),  con- 
sumption or  energy  requirements  in  any  other  terms 
— such  as  kilowatt  hours  supplied  to  or  required  by  a 
centrifugal  chiller.  Btu’s  of  chilled  water  required, 
gallons  of  oil  or  cubic  feet  of  gas  supplied  to  a boiler, 
or  gallons  of  hot  water  used — must  be  converted  to 
thermal  energy  in  the  form  of  hot  water  (or  hot  air  if 
a hot  air  system  is  used).  This  usually  requires  taking 
into  account  the  efficiency  of  conversion  equipment 
as  well  as  the  appropriate  conversion  factors.  Deter- 
mination of  Ql  for  the  three  types  of  systems  is 
described  below. 

Solar  Domestic  Hot  Water  System.  To  use  the 

universal  curve  for  domestic  hot  water  systems,  the 
annual  energy  needed  for  heating  the  hot  water  is  re- 
quired. The  best  source  of  information  is  metered 


As  the  example  in  the  previous  section  shows,  two 
sets  of  input  information  in  addition  to  collector  area 
are  required  to  use  the  universal  curve:  the  thermal 
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data  from  the  building  in  question  or  from  an  identi- 
cal or  similar  building  at  the  same  location.  Metered 
consumption  data  for  either  energy  or  fuel  input  to 
the  hot  water  heater  or  gallons  of  hot  water  con- 
sumed can  be  used.  When  fuel  consumption  is 
known,  corrections  for  boiler  efficiency  must  be 
taken  into  account  as  discussed  above.  When  actual 
consumption  data  cannot  be  obtained,  data  can  be 
derived  from  various  handbooks  or  from  local  or 
national  plumbing  codes.  However,  care  must  be 
exercised  in  using  data  in  the  torm  ot  per  capita  con- 
sumption, as  they  are  usually  developed  for  design 
purposes  and  thus  are  often  more  closely  related  to 
peak  demand  than  to  average  hot  water  consump- 
tion. Whichever  method  is  used,  if  annual  hot  water 
demand  is  determined,  'he  annual  energy  consump- 
tion is  obtained  from  Eq  2: 

Ql  = DwCp  (Tou,  — Tjn)  IF.q  2| 

where  Qp  = annual  energy  consumption  in  Btu(J) 
D = annual  hot  water  demand  in  gal(f) 
w = density  of  water  = 8.33  lb  gal 
(1.00  kg/£) 

C„  = specific  heat  of  water  = 1 Btu/lbm 
°F  (4.1868  kJ /kg  °K) 

Tout  = temperature  of  hot  water  supply  in  °F 
<°C) 

Tjn  = temperature  of  supply  water  in  °F  (°C) 

Solar  Heating  Systems.  Since  the  length  of  the 
heating  season  varies  significantly  from  site  to  site, 
the  universal  curve  for  solar  heating  systems  (liquid 
or  air)  must  be  used  on  a monthly  basis.  This  in  turn 
requires  that  the  monthly  energy  demand  for  heating 
be  obtained. 

Again,  the  best  source  of  data  for  estimating  the 
amount  of  heat  energy  required  for  a particular 
building  or  set  of  buildings  is  measured  data  for  the 
building  in  question  or  for  similar  or  identical  build- 
ings at  the  same  location.  These  data  may  only  be 
available  in  terms  of  the  amount  of  fuel  oil  used  per 
month  or  the  amount  of  gas  consumed  in  heating  the 
particular  building  or  buildings.  If  so,  the  fuel  con- 
sumption data  must  be  adjusted  to  account  for  the 
efficiency  of  the  boiler  or  furnace  supplying  the  heat 
to  the  space.  Typically,  solar  heating  systems  include 
both  heating  the  building  and  the  domestic  hot  water 
for  the  building.  This  can  be  accounted  for  in  the 
energy  requirement  by  adding  the  monthly  domestic 
hot  water  load  to  the  heating  load. 

If  measured  data  for  the  building  being  con- 
sidered are  not  available,  the  building's  heating 


energy  demand  can  be  estimated  using  one  of  several 
methods.  One  such  method,  described  in  a CKRL 
report  entitled  Energy  Utilization  Method  for  Pre- 
dicting Pudding  Energy  Use.3  provides  a monthly 
energy  utilization  index  in  Btu’s  per  square  foot  per 
month  for  several  typical  Army  buildings  and  a 
method  for  adjusting  this  estimate  based  on  the 
differences  between  the  building  being  considered 
and  the  typical  buildings  for  which  the  energy  use 
indices  were  developed.  The  energy  use  index 
method  can  also  be  used  in  conjunction  with  meas- 
ured data  for  a building  similar  to  the  one  being  con- 
sidered by  adjusting  the  data  for  the  similar  building 
based  on  the  ratio  of  the  two  buildings’  energy 
indices. 

The  Carrier  Rational  Energy  Analysis  Procedure 
(REAP)4  is  another  method  for  calculating  building 
heating  requirements.  This  method  is  a modified  bin 
method5  which  requires  the  calculation  of  heating 
loads  at  several  different  design  points  of  the  system. 
These  loads  combined  with  frequency  of  occurrence 
of  weather  conditions,  which  can  be  obtained  from 
Air  Force  Manual  88-8, 6 provide  an  estimate  of  the 
monthly  energy  requirements  for  heating. 

A fourth  method  requires  application  of  the  Build- 
ing Loads  Analysis  and  System  Thermodynamics 
(BLAST)  computer  simulation  program  to  predict 
the  monthly  heat  energy  demands.  Use  of  this  pro- 
gram, which  is  described  in  the  program  user's 
manual.7  requires  an  hourly  weather  data  tape. 

Monthly  energy  consumption  for  domestic  hot 
water  heating,  which  is  generally  included  in  a solar 
heating  system,  can  be  determined  by  applying  the 
methods  described  in  the  preceding  section  on  a 
monthly  basis.  It  should  be  noted  that  measured 
data  often  have  domestic  hot  water  and  space  heat- 
ing combined. 


’L.  M.  Windingland  and  D.  C.  Hittle.  Energy  Utilization 
Index.  Method  for  Predicting  Building  Energy  Use.  Draft  Techni- 
cal Report  (CERL.  1977). 

4 Rational  Energy  Analysis  Procedure  ( REAP ) (Carrier  Air  Con- 
ditioning Company,  undated). 

5 ASHRAE  Handbook  of  Fundamentals  (American  Society 
of  Heating.  Refrigerating  and  Air  Conditioning  Engineers 
| ASHRAE],  1972). 

4 Engineering  Weather  Data.  AEM  88-8.  Chapter  6 (Depart- 
ment of  the  Air  Force.  1967). 

7Building  Loads  Analysis  and  System  Thermodynamics 
{BLAST)  Program  Users  Manual.  Draft  Technical  Report 
(CERL.  1977). 
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Design  day  calculations  alone  cannot  be  used  to 
obtain  heating  loads,  since  they  provide  peak  energy 
demands  rather  than  the  average  demand  required 
for  estimating  monthly  energy  consumption. 

Solar  Heating  and  Cooling  Systems.  I he  universal 
curse  for  solar  heating  and  cooling  systems  is  an 
annual  curve,  since  there  is  a year-round  demand  for 
energy.  Thus,  the  annual  energy  is  needed.  The 
annual  energy  required  is  the  sum  of  the  heating  and 
cooling  energy  requirements  plus  the  domestic  hot 
water  heating  energy  requirement  if  it  is  to  be  in- 
cluded in  the  system. 

The  annual  energy  required  for  heating  and  cool- 
ing is  again  best  determined  from  measured  data  for 
the  building  in  question  or  a similar  or  identical 
building  at  the  same  location.  In  determining  the 
energy  required  for  cooling,  attention  must  be  given 
to  the  differences  in  the  coefficients  of  performance 
(COP)  of  the  various  chillers  involved.  For  example, 
if  measuring  the  consumption  of  electricity  for  a con- 
ventional chiller  in  the  building  indicates  that 
1 .(XX), 000  kWh  are  used  annually,  then  this  number 
must  first  be  converted  to  Btu’s  by  multiplying  by 
3412,  giving  3.41  x 109  But  (or  converted  to  kilo- 
joules by  multiplying  by  3600.  giving  3.60  x 109  kJ). 
and  then  multiplied  by  the  centrifugal  chiller’s  mean 
COP.  If  the  mean  COP  is  4,  13.6  x 109  Btu/yr  (14.3  x 
109  kJ/yr)  are  required  to  be  delivered  by  the  chiller 
to  the  cooling  system.  This  chilled  water  demand, 
when  divided  by  the  COP  of  a solar  energy  system's 
absorption  chiller  (say  .65)  indicates  that  21.0  x 109 
Btu  (22.2  x 109  kJ)  must  be  supplied  to  the  absorp- 
tion chiller  to  meet  the  annual  energy  required  for 
cooling.  Since  heating  and  cooling  systems  may  also 
include  domestic  hot  water  heating,  total  demand 
for  heat  energy  for  heating  domestic  hot  water  can  be 
summed  with  the  total  demand  for  heating  and 
added  to  the  annual  demand  for  cooling  energy  in 
order  to  determine  the  annual  building  energy  re- 
quirement. Ol- 

If  measured  consumption  data  are  not  available, 
the  energy  use  index  method,  REAP  method,  or 
BLAST  simulation  program  described  in  the  pre- 
vious section  can  be  used  to  determine  annual  heat- 
ing and  cooling  loads. 

Solar  Radiation  Data 

The  tilt  angle  of  the  collector  plate  greatly  affects 
the  amount  of  solar  radiation  striking  the  surface  of 
the  plate.  Most  measured  and  reported  data  are  for 


solar  radiation  striking  a horizontal  surface.  This 
section  describes  where  to  obtain  this  inl<  nation 
and  how  to  correct  it  for  tilted  collector  arrays 

Average  solar  radiation  data  for  many  sites 
around  the  country  are  published  in  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA) 
Climatic  Atlas  of  the  United  States.  Annual  and 
monthly  solar  radiation  maps  from  this  document 
are  reproduced  in  Appendix  A.  Since  the  radiation 
values  on  these  maps  are  given  in  daily  means,  the 
annual  and  monthly  values  are  obtained  by  multiply- 
ing by  the  number  of  days  in  the  year  or  month, 
depending  on  the  map  being  used. 

The  Climatic  Atlas  of  the  United  Slates  also  con- 
tains tabulated  summaries  ot  radiation  data  for 
specific  sites.  More  detailed  summaries  are  fre- 
quently available  from  NOAA  or  from  local  weather 
services  or  state  agencies.  Data  for  the  particular  site 
in  question  should  be  obtained,  if  possible,  since 
incident  solar  radiation  frequently  varies  substan- 
tially over  relatively  small  geographical  distances. 
This  is  particularly  true  in  coastal  regions  and  in 
regions  near  mountains,  where  local  climatological 
variations  are  severe. 

Once  the  collector  tilt  angle  is  known,  the  radia- 
tion flux  density  on  the  filteci  collector  surface  can  be 
determined  from  horizontal  radiation  data.  Fortu- 
nately. there  is  an  optimum  collector  tilt  angle 
(measured  from  horizontal)  fpr  each  type  of  system. 
For  heating  and  cotiling  systems,  this  optimum  tilt 
angle,  0C.  is  roughly  the  location  latitude  minus  10 
degrees;  for  heating-only  systems,  it  is  roughly  the 
latitude  plus  10  degrees;  and  for  domestic  hot  water 
systems,  it  is  roughly  equal  to  the  latitude.  These 
angles  provide  for  collection  of  the  greatest  amount 
of  olar  energy  in  each  application.  Variations  of  ±5 
degrees  affect  performance  only  slightly.  In  all  cases 
the  optimum  azimuth  angle  is  due  south  and  again 
slight  deviations  from  due  south  (±10  degrees)  do 
not  significantly  reduce  system  performance. 

Once  the  optimum  tilt  angle  has  been  determined, 
the  annual  radiation  tlux  density,  Qc.  on  the  opti- 
mally tilted  surface  can  be  estimated  from  Eq  3: 

Qc  = KHjj  [Eq  3] 

where  K = a correction  factor  depending  on  the 
collector  tilt 

= the  annual  or  monthly  radiation  flux 
density  on  a horizontal  surface  (Btu/ 
sq  ft  or  J/m2). 
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For  solar  systems  for  which  annual  solar  radiation 
data  are  used  (i.e. . heating  and  cooling  or  domestic 
hot  water),  the  correction  factor,  K,  is  given  approxi- 
mately by 


K = 


cos(i?L  -7  — t?c) 
cos  (i?l  ~7) 


[ Ecj  4 1 


where  i?l  = the  latitude  in  degrees 

dc  = the  optimum  angle  from  the  horizontal 
in  degrees. 

Eq  4 is  an  empirically  derived  equation  based  on 
the  results  of  simulation  studies  and  is  valid  only  for 
near-optimum  collector  tilt  angles  and  annual  solar 
radiation  data. 

For  solar  heating  applications,-,  where  monthly 
radiation  figures  must  be  used,  the  correction  factor 
may  not  be  expressed  in  closed  form.  In  this  case.  K 
must  be  determined  from  Table  1.  which  gives 
monthly  values  of  this  correction  factor  for  four 
different  latitudes.  Table  1 is  valid  only  for  optimum 
collector  tilt  angles. 

Table  1 

Monthly  Correction  Factors  IK)  for  Heating-Only  Systems 
(for  collector  tilt  equal  to  latitude  plus  10  degree*) 


Month 

30° 

Latitude 
35°  40° 

45° 

Jan 

1.64 

1.84 

2.12 

2.49 

Feb 

1.42 

1.55 

1.71 

1.93 

Mar 

1.18 

1.25 

1.33 

1.44 

Apr 

.97 

1.00 

1.03 

1.07 

Mav 

.84 

.84 

.85 

.87 

Jun 

.78 

.78 

.78 

.78 

Jut 

.80 

.80 

.80 

.81 

Aug 

.90 

.91 

.93 

.96 

Sep 

1.07 

III 

1. 17 

1.24 

Oct 

1.29 

1.40 

1.51 

1.66 

Nov 

1.54 

1.70 

1.93 

2.23 

Dec 

1.70 

1.93 

2.24 

2.64 

Recalling  the  example  on  page  8,  the  annual  radi- 
ation for  Topeka,  KS  was  found  as  follows:  the 
annual  radiation  level  on  a horizontal  surface  ob- 
tained from  Figure  A1  is  380  langleys/day.  The  lati- 
tude. i?l  (also  *rom  Figure  A 1 ).  is  39.  Using  Eqs  3 
and  4.  Qc  can  be  found. 


COsftJl 


<V 


_ cos  (39  -7-39)  /380  langleysA 
e cos(39  — 7)  V day  / 

/365day\  / 3.69  Btu  \ 

\ yr  / \ sq  ft-langley  ) 


°c  = »<H,=  cos(,l-_7)  ", 


Qc  = 6.0  x 10s  Btu/sq  ft/yr 
(6.8  x 106kJ/ml/yr) 

As  an  example  of  applying  the  monthly  correction 
factor  used  for  heating-only  systems,  find  the  solar 
radiation  for  January  in  Topeka.  KS.  From  Table  1, 
K is  estimated  by  extrapolation  between  35  degrees 
and  40  degrees  latitude  for  January.  Thus 

K-|84+(?#^)#'-)x4' 

K =2.06 

From  Figure  A2,  the  mean  daily  horizontal  radiation 
is  190  langleys.  Hence, 

H _ / 190  langleys  \ / 31  days\  / 3.69  Btu  \ 

V day  / \ mo  / V sq  ft-langley  / 

= 2.2  x 104  Btu/sq  ft/mo  (2.5  x 10*  kJ/m2/mo) 

Thus, 

Qc  = KHfl  = 4.5x|04  Btu/sq  ft /mo 
(5. 1 x 10s  kJ/m2/mo) 


ESTIMATING  SOLAR  SYSTEM 
PERFORMANCE 


General 

In  determining  optimal  collector  area  and  making 
the  energy  and  life-cycle  cost  analysis,  a life-cycle 
cost  comparison  between  a conventional  system  and 
the  solar  energy  system  for  various  collector  areas  is 
made.  A critical  step  in  this  analysis  is  determination 
of  the  solar  system  performance  to  find  out  how- 
much  of  the  energy  requirement  is  met  by  solar 
energy  and  how  much  is  met  by  auxiliary  heat  for  a 
given  size  collector  array.  In  order  to  make  this 
determination,  the  universal  curves  are  used. 

Analyzing  a given  solar  system  performance  re- 
quires the  following  steps: 
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1.  Determine  the  annual  or  monthly  energy  re- 
quirement. as  described  on  pages  9-11. 

2.  Determine  the  annual  or  monthly  solar  radia- 
tion on  the  optimally  tilted  collector  array,  Qe,  lor 
the  location  in  question  using  the  method  described 
on  pages  11-12. 

3.  Determine  the  solar  system  performance  using 
the  universal  curve.  Two  methods  can  be  used  to 
accomplish  this: 

a.  If  a collector  area,  Ac,  is  assumed,  then  the 
annual  or  monthly  solar  system  performance  param- 
eter, Ps,  can  be  found  from  Eq  1: 

0CAC 

ps=75T  lEq1' 

From  this  value  and  the  universal  curve,  the  annual 
or  monthly  fraction  of  energy  met  by  the  solar 
system,  p,  can  be  determined.  Using  this  value  in 
Eqs  5 and  6 yields  the  energy  supplied  by  auxiliary 
heat,  OlA'  an<*  t*le  ener8y  supplied  by  the  solar 
system , °LS 

QLA=«-e)QL  (Eq  6] 

If  the  analysis  is  for  heating-only  systems,  the 
monthly  QlS  and  (?LA  are  obtained  from  Eqs  5 and 
6 and  must  be  summed  to  get  the  annual  values  re- 
quired for  determining  optimal  collector  area  and 
making  the  energy  and  life-cycle  cost  analysis  as 
described  in  Chapter  4. 

b.  The  universal  curve  can  also  indicate  what 
collector  area  is  required  to  meet  a given  fraction,  q. 
of  the  annual  or  monthly  energy  requirement.  This 
information  is  useful  in  establishing  a starting  point 
for  a solar  system  analysis.  The  collector  area  can  be 
determined  by  using  q to  find  Ps  from  the  universal 
curve.  For  Ps.  the  collector  area  can  be  found  from: 


Again,  annual  or  monthly  Qls  and  QlA  can  be 
found  from  Eqs  5 and  6. 

The  next  sections  present  the  universal  curves  for 
the  various  systems  and  examples  of  how  to  use 
them. 


Solar  Domestic  Hot  Water  Systems 

Development  of  the  universal  curve  for  domestic 
hot  water  heating  (Figure  1)  assumed  a uniform 
annual  hot  water  demand.  If  the  system  being 
analy/.ed  has  a highly  seasonal  demand,  this  curve 
should  not  be  used.  The  curve  also  was  derived  for  a 
system  which  has  an  inlet  water  temperature  of  55°F 
(I3°C)  and  a supply  temperature  of  140°F  (60°C), 
which  are  typical  of  hot  water  systems.  The  curve  is 
valid,  however,  for  small  variations  (±20°F  (±1 1 °CJ) 
about  these  values. 

The  following  example  illustrates  application  of 
the  procedure  outlined  above  for  solar  domestic  hot 
water  heating. 

Example  1: 

Assume  a barracks  at  Fort  Hood,  TX  consumes 
30  gal  (114,0  of  hot  w ater  per  day  per  person  and  the 
average  occupancy  for  the  year  is  200  people.  The 
collector  area  which  will  meet  50  percent  of  the  load 
if  the  city  supply  water  is  55°F  (13°C)  and  the  hot 
water  system  supplies  140°F  (60°C)  water  can  be 
determined  in  the  following  manner: 

1.  First  determine  Ql-  the  annual  energy  re- 
quired to  heat  the  water.  Eq  2 gives: 

0L  = DwCp(Tout-Tin) 

°l  - (“far,)  1200 Cv”) 

(8-fr)  ( w)  ll«"F'55"F» 

QL  = 1.6*  109  Btu/yr  (1 .7  x 109kJ/yr) 

2.  Estimate  Qc,  the  annual  solar  radiation  on  the 
array.  From  Eqs  3 and  4, 

cos  (0L  - 7 - Cc) 

°c=  cos((7l-7)  H0 

From  Figure  Al,  the  mean  daily  horizontal  solar 
radiation  at  this  location  is  445  langleys/day.  Thus, 

_ / 445  langleys\  / 365days\  / 3,69  Btu  \ 
® \ day  /\  yr  / \sq  ft-langley  / 

= 6.0  x 10s  Btu/sq  ft /yr  (6.8  x 10*  kJ/mVyr) 
Also  from  Figure  Al.  i?l  is  33  degrees  (latitude  of 
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Fort  Hood.  TX).  Since  this  is  a domestic  hot  water 
system,  the  optimum  tilt  angle.  t)t..  for  hot  water  is 
equal  to  i>l  (33  degrees).  Thus. 

Qc  =6.6  x 10s  Btu  ^sq  ft/yr  (7.5  x 10*  Kj/ml/yr) 

3.  Use  the  universal  curve  to  estimate  the  solar 
performance.  To  meet  50  percent  of  the  load,  q is  .5. 
From  Figure  1 (for  q equal  to  .5)  ^ is  .85.  The  re- 
quired collector  area  can  now  be  found  from  Eq  I. 

A =PsQl  (.85)11.6  x 10*  Btu /yr) 

C Qc  (6.6  x 10*  Btu  sq  ft-yr) 

= 2 x 101  sq  ft  (186  ml) 

The  area  required  to  meet  50  percent  of  the  annual 
domestic  hot  water  load  is  2000  sq  ft  (186  mJ). 


Solar  Heating  and  Cooling  Systems 

Figure  2 shows  the  universal  curve  for  solar  heat- 
ing and  cooling  systems.  Analysis  of  the  system’s 
performance  proceeds  exactly  as  described  in  the 
previous  section  for  domestic  hot  water  systems.  The 
following  example  illustrates  the  application  of  the 
solar  heating  and  cooling  universal  curve. 


Example  2: 

The  following  example  is  for  an  administrative 
building  housing  300  people  in  Champaign,  IL  (lati- 
tude 42  degrees).  The  energy  consumption  data  were 
taken  from  monthly  utility  billings.  The  steps  in  the 
analysis  are  as  follows: 

1 . First  determine  Q^.  the  annual  energy  require- 
ment for  heating  and  cooling. 

a.  For  heating,  power  records  indicate  that  for  the 
months  of  October  through  June,  144  x 10*  Btu 
(152  x 10*  kJ)  of  gas  were  consumed  for  heating  and 
domestic  hot  water.  When  multiplied  by  a boiler 
efficiency  of  .75,  the  building  heating  load  is  found 
to  be  108  x 10*  Btu  (114  x 10*  kJ). 

b.  For  cooling,  utility  records  show  that  for  the 
months  of  April  through  November,  81  x 104  kWh 
(28  x 10*  Btu  |29.5  x 10*  kJ))  were  consumed  by  the 
centrifugal  chiller.  The  chiller  load  was  found  by 
subtracting  the  average  winter  electrical  load  from 
the  average  summer  electrical  load  to  account  for 
other,  noncooling  loads.  Thus,  chiller  electrical 
load,  when  multiplied  by  the  COP  of  the  centrifugal 
chiller  (found  from  manufacturer's  literature  to  be 
4),  gives  a cooling  load  for  the  building  of  112  x 10* 
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ANNUAL  INCIDENT  RADIATION  / ANNUAL  ENERGY  REQUIREMENTS,  Ps 
Figure  2.  Universal  curve  for  solar  heating  and  cooling  systems. 


Btu  (118  x 10*  kJ).  Dividing  by  the  COP  of  the 
absorption  chiller  (.65)  indicates  that  172  x l()*  Btu 
(181  x 10*  kJ)  must  be  supplied  to  the  absorption 
chiller  to  meet  the  annua!  cooling  energy  demand. 

Hence. 

Ql  = 108  x 10*  Btu/yr  + 172  x 10*  Btu/yr 
0L  = 280  x 10*  Btu/yr  (295  x 10*  kJ/yr) 

2.  Estimate  the  annual  solar  radiation.  Qc,  on  the 
optimally  tilted  collector  array.  From  Eqs  3 and  4 

cos(i?i  — 7-i?J 

o _ u 

c cos(tJ^  — 7)  0 

From  Figure  Al,  the  mean  daily  annual  solar  radia- 
tion for  Champaign,  IL  is  355  langleys/day.  Thus 
Hjj  is 

^ _ 355  langleys  365  days  3.69  Btu 

» day  yr  sq  ft-langley 

= 4.8  x |05  Btu/sq  ft/yr  (5.5  x 10*  kJ/m2/yr) 

since  i7^  =42  degrees  (the  latitude  of  Champaign) 

and  i?c  = — 10  degrees  = 32  degrees 

(optimum  collector  tilt  angle  for  heating 
and  cooling) 

Thus,  Qc  = 5.6  x 105  Btu/sq  ft/yr 
(6.4  x 10‘kJ/mVyr) 

3.  Determine  the  solar  system  performance  from 
the  universal  curve  for  heating  and  cooling  (Figure  2) 
assuming  a collector  area  of  50,000  sq  ft  (4645  m2). 
From  Eq  1,  Ps  is  given  by 

_ _ ^c^c  _ (5  x 104sq  ft)(5.6  x 10s  Btu/sq  ft/yr) 

= QL“  = (280x16*  Btu/yr) 

Ps  = 1.0 

For  Ps  = 1,  Figure  2 shows  q = .26. 

Hence,  roughly  26  percent  of  the  building’s  annual 
energy  requirement  for  heating  and  cooling  can  be 
supplied  by  the  sun  using  50,000  sq  ft  (4645  m2)  of 
collectors  for  an  annual  energy  savings  of  78.4  x 10* 
Btu/yr  (82.7  x 10«  kJ/yr). 


Solar  Heating  Systems 

Figure  3 shows  the  universal  curve  for  solar  heat- 
ing systems.  Methods  for  using  (he  curve  are  similar 
to  those  for  the  two  previous  systems,  except  the 
analysis  is  performed  month  by  month. 

Since  the  domestic  hot  water  load  for  a heating 
system  is  small  in  comparison  to  the  space  heating 
requirements,  it  can  be  assumed  that  the  total  hot 
water  load  is  met  by  solar  energy  during  the  non- 
heating season.  This  somewhat  simplifies  the  calcu- 
lations, in  that  the  monthly  analysis  needs  to  be  per- 
formed only  for  those  months  with  significant  space 
heating  loads.  The  hot  water  energy  requirements 
for  the  other  months  can  be  summed  and  added 
directly  to  Ols-  The  following  example  shows  how 
this  monthly  analysis  is  accomplished  for  heating- 
only  systems. 

Example  3: 

The  administrative  building  of  Example  2 was 
used  in  this  example  to  determine  what  percent  of 
the  heating  load  could  be  met  by  the  50,000  sq  ft 
(4645  m2)  of  collectors.  Again  the  energy  consump- 
tion data  were  taken  from  monthly  utilities  bills.  The 
steps  in  the  analysis  are  as  follows: 

1.  Determine  the  monthly  loads,  Ol-  f°r  the 
heating  season. 

Table  2 presents  the  building's  natural  gas  con- 
sumption indicated  by  the  monthly  utility  records. 
The  values  of  Ql  in  the  table  have  been  corrected  for 
a boiler  efficiency  assumed  as  .75.  Note  that  the 
108  x 10*  Btu/yr  (114  x 10*  kJ/yr)  total  load  is  the 
same  as  found  in  Example  2. 

2.  Estimate  Qc.  the  annual  solar  radiation  on  the 
optimally  tilted  collector  array.  From  Eq  3 

Qc  = KHfl 

where  H^  and  Qc  = monthly  radiation  values. 

K = the  correction  factor  from  Table  I. 

Hfl  is  obtained  from  the  figures  in  Appendix  A. 
Monthly  values  for  Qc  are  listed  in  Table  3. 

3.  Determine  the  solar  system  performance  from 
the  universal  curve  for  heating  only  (Figure  3), 
assuming  a collector  area,  Ac,  of  50,000  sq  ft  (4645 


Table  2 Table  3 


Natural  Gas  Consumption  for  Example  Administration  Building  Monthly  Values  for  Qc 


(Heating) 

Ql(  Hot  Water) 

Total 

no 

Qc 

Month 

Btu 

x 10* 

Btox  10* 

Btu/ mo  x 10* 

Month 

Langley/day 

Btu/sq  ft/mo 

K 

Btu/sq  ft/mo 

(Id* 

10«l 

(Id*  10*) 

(kj/tno) 

(Id/m'/mo) 

(U/mVmo) 

Oct 

2.0 

(2.1) 

.6  (.6) 

2.6 

(2.7) 

Nov 

7.8 

(8.2) 

.6  (.6) 

8.4 

(8.9) 

Oct 

275 

3.1  x 10* 

1.5 

4.7x  10* 

Dec 

12.3 

(13.0) 

.6  (.6) 

12.9 

(13.6) 

<3.5x  10») 

(5.3  x 10s) 

Jan 

19.9 

(21.0) 

.6  (.6) 

20.5 

(21.6) 

Nov 

175 

1.9x  10* 

2.0 

3.8  x 10* 

Feb 

18.5 

(19.5) 

.6  (.6) 

19. 1 

(20.2) 

(2.2x  10s) 

(4.3x  10s) 

Mar 

16.3 

(17.2) 

.6  (.6) 

16.9 

(17.8) 

Apr 

12.9 

(13.6) 

.6  (.6) 

13.5 

(14.2) 

Dec 

135 

1.5x  10* 

2.3 

3.5*  10* 

May 

8.8 

(9.3) 

.6  (.6) 

9.4 

(9.9) 

(l.7x  10s) 

(4. Ox  10s) 

Jun 

4.0 

(4.2) 

.6  (.6) 

4.6 

(4.9) 

Jan 

155 

!.8x  10* 

2.2 

4.0  x 10* 

(2. Ox  10s) 

(4.5  x 10s) 

102. 5«  10*  Btu/yr 

I07.9x  10*  Btu/yr 

2.5  x 10* 

1 108. 1 x 10*  kj/yr) 

( 1 1 3.8 x 10*  kj/yr) 

Feb 

240 

1.8 

4.5x  10* 

(2.8  x 10*) 

(5.1  x 10s) 

m2).  From  Eq  1 , Ps  for  each  month  can  be  calculated 

Mar 

330 

3.8  x 10* 
<4.3  x 10s 

1.3 

4.9x  10* 
(5.6  x 10s) 

trom  the  monthly  Qi  and  Qc  presented  in  Tables  2 

Apr 

400 

4.4  x I04 

1.0 

and  3. 

4.4  x 10* 

(S.Ox  10s) 

(S.Ox  10s) 

From  these  values  of  P.-.  the  monthly  values  of  o 

May 

510 

5.8  x 10* 

.85 

4.9x  10* 

can  then 

be 

obtained 

from  the  universal 

curve. 

(6.6  x 10s) 

(5.6  x 10s) 

Finally,  froirfEq  5,  the  load  met  by  the  solar  system, 
Qi  s-  can  be  found  for  each  month  (Table  4). 

lun 

550 

6.1  x 10* 
(6.9  x 10s) 

.78 

4.8x  10* 
(5.5  x 10s) 
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Table  4 


Monthly  Value*  of  e and 


Month 

<?c 

Btu/mo  Btu/aq  ft/mo 
(ld/mo)  IkJ/mVmo) 

P. 

0 

qls 

Btu 

(Id) 

Oct 

2.6  » 10* 
(2.7  x 10*) 

4.7x  10* 
(5.3  x 10*) 

9.04 

.99 

2.6  x 10* 
(2.7  x 10*) 

Nov 

8.4  x 10* 
(8.9x  10*) 

3.8x  10* 
(4.3  x 10s) 

2.26 

.59 

S.Ox  10* 
(S.2x  10*) 

Dec 

12.9x  10* 
<13.6  x I01) 

3.5x  10* 
(4.0  x 10s) 

1.36 

.39 

S.Ox  10* 
(5.2  x 10*) 

Jan 

20.5x  |0» 
(2l.6xl0«) 

4. Ox  10* 
<4.Sx  10s) 

.98 

.29 

5.9  x 10* 
(6.3  x 10*) 

Feb 

19.1  x 10* 
(20.2  x 10*) 

4.5x  10* 
(5. 1 x 10s) 

1.18 

.35 

6.7x  10* 
(7.0  x 10*) 

Mar 

I6.9x  10* 
(17.8x  10*) 

4.9x  10* 
(5.6  x 10s) 

1.45 

.41 

6.9x  10* 
(7.3x  10*) 

Apr 

13.5  x 10* 
(14.2  * 10*) 

4.4  x 10* 
(5.0  x 10s) 

1.63 

.45 

6.1  x 10* 
<6.5x  10*) 

May 

9.4  x 10* 
(9.9  x 10*) 

4.9x  10* 
(5.6  x 10s) 

2.61 

.65 

6. 1 x 10* 
(b.5x  10*) 

Jun 

4.6  x 10* 
(4.9  x 10*) 

4.8x  10* 
(5.5  x 10s) 

5.22 

.97 

4.5x  10* 
(4.7  x JO*) 

Total 

°LS 

48.8  x 10*  Btu 

(51.6X  * 0*kJ) 

The  figures  in  Table  4 show  that  48.8  x 10*  Btu 
(51.6  x 10*  kJ)  of  heating  is  supplied  by  the  solar 
system  during  October  through  June.  The  total  heat- 
ing Of  f°r  October  through  June  (from  part  1 of  this 
example)  is  108  x 10*  Btu  (1 14  x 10*  kJ).  Hence,  the 
solar  system  provides  45  percent  of  the  annual  heat- 
ing requirement. 

For  the  3 months  not  considered  above  it  is 
assumed  that  all  domestic  hot  water  heating  is  per- 
formed by  the  solar  energy  system.  Thus, 

Qls  = 48.8  x 10*  Btu  + 3(.6  x 10*  Btu) 

= 50.6  x 10*  Btu  (63.5  x 10*  kJ) 

and 

Ql  = 107.9  x 10*  Btu  + 3(.6  x 10*  Btu) 

= 109.7  x 10'  Btu  (1 15.7 x 10' kJ) 
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Thus,  the  50,000  sq  ft  (4645  ml)  of  collectors  will 
meet  46  percent  of  the  heating  and  hot  water  load  for 
a total  energy  savings  of  50.6  x 10*  Btu  (53.4  x 10* kJ) 
compared  with  Example  2,  where  the  same  50,000 
sq  ft  (4645  m1)  of  collectors  met  only  26  percent  of 
the  heating  and  cooling  load  but  had  a total  energy 
savings  of  78.4  x 10*  Btu  (82.7  x 10'  kJ). 

Corrections  for  Different  Collector  Types 

As  previously  described,  the  universal  curves  were 
developed  for  the  reference  solar  collector.  For 
collectors  with  different  absorber  plate  characteris- 
tics and  different  numbers  of  covers,  Tables  5 and  6 
can  be  used  to  correct  the  results  of  the  universal 
curve.  These  tables  give  values  for  a collector  area 
multiplication  factor  which  can  be  applied  when  a 
particular  collector  other  than  the  reference  collector 
is  used.  In  these  tables,  a is  the  absorptivity  and  t is 
the  emissivity  of  the  absorber  plate.  The  transmit- 


TableS 

Collector  Are*  Multiplying  Factor*  for  Different  Collector 
Design!  for  Heating  and  Cooling  Systems 


a 

0.96 

0.94 

0.90 

c 

0.96 

0.30 

0.10 

N 

i 

1.55 

1.09 

1 

2 

1.09 

0.97 

0.93 

a = Absorptivity 
t — Emissivity 
N Number  of  glass  covers 


Table  6 

Collector  Area  Multiplying  Factor*  for  Different  Collector 
Design*  for  Domestic  Hot  Water  Heating 


and  Heating-Only  Syatem* 

a 

.96 

.90 

c 

.96 

.10 

N 



1 

1.26  1 

2 

.96 

.93 

Hence,  the  annual  q is 
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a = Absorptivity 

c = Emissivity 

N - Number  of  glass  covers 


lance  ot  each  cover  was  assumed  to  b 0.9,  so  that  for 
two  covers,  the  transmittance  equals  0.81.  These 
tables  are  useful  because  they  permit  application  of 
the  performance  curves  for  the  reference  flat-plate 
collector  (a  single-cover  collector  with  an  absorp- 
tivity of  0.9  and  emissivity  of  0.10)  to  other  flat -plate 
collectors.  For  example,  in  order  to  achieve  equiva- 
lent performance  for  heating  and  cooling  from  a two- 
cover.  nonselective  collector,  one  would  multiply  the 
appropriate  collector  area  for  the  reference  collector 
by  1.09  or  increase  the  collector  area  by  9 percent. 
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DETERMINING  THE  OPTIMAL 
COLLECTOR  AREA 


This  chapter  describes  a method  for  defining  the 
optimal  collector  area  which  is  based  on  a life-cycle 
cost  comparison  of  the  solar  energy  system  and  a 
conventional  energy  system.  The  procedure  is 
adapted  from  that  suggested  by  Butz,  et  al.' 

The  method  relies  on  a simple  graphical  compari- 
son of  the  life-cycle  cost  difference  between  a solar 
energy  system  and  a conventional  energy  system, 
where  the  solar  system  performance  is  determined 
from  the  universal  curve  as  described  in  Chapter  3. 
The  Office  of  the  Chief  of  Engineers  (OCE)  life-cycle 
cost  instructions’ should  be  used.  Since  the  purpose 
of  this  optimization  method  is  to  determine  the  eco- 
nomically optimum  collector  array,  life-cycle  cost 
differences  are  examined  rather  than  total  system 
life-cycle  costs.  Normal  procedures  and  instructions 
for  determining  overall  building  and  energy  system 
capital  and  life-cycle  costs  of  various  options  can  be 
followed  once  the  optimum  system  is  determined. 

The  first  step  in  determining  the  optimum 
collector  area  is  to  estimate  the  life-cycle  cost  of  the 
solar  energy  system  components.*  Because  the 
method  being  described  is  comparative,  only  the  cost 
of  components  that  are  not  normally  part  of  the  con- 
ventional heating  and  cooling  system  should  be  con- 
sidered. For  example,  the  cost  of  the  building's  air- 
handling system  would  not  be  considered,  but  the 


•L  . W.  But/,  et  at.,  Use  of  Solar  Energy  for  Residential  Healing 
and  Cooling.  Masters  Thesis  in  Mechanical  Engineering  (Univer- 
sity of  Wisconsin.  1973). 

'OCE  Life  Cycle  Costing  Instructions  (Department  of  the 
Army.  May  1971). 

•This  report  does  not  include  cost  data,  since  solar  component 
costs  are  changing  rapidly.  Manufacturers  are  the  best  source  of 
data. 


difference  between  the  cost  of  a more  expensive 
absorption  chiller  and  a less  expensive  centrifugal 
chiller  should  be  charged  to  the  solar  energy  system. 

Certain  cost  elements  for  solar  energy  systems  vary 
according  to  the  size  of  the  solar  energy  system, 
while  others  are  relatively  fixed  regardless  of  collec- 
tor area  or  tank  size.  Collector  and  storage  tank 
costs  are  obvious  examples  of  costs  which  are 
dependent  on  collector  area.  (The  tank  volume  is 
assumed  to  be  16  lb  of  water  per  square  foot  of  col- 
lector [78kg/m,|.)  Other  examples  include  heat  ex- 
changer costs  and  certain  pump  and  piping  costs. 
Additional  control  system  costs  associated  with  the 
solar  energy  system  are  an  example  of  cost  differ- 
ences which  are  largely  independent  of  collector 
area.  The  cost  difference  associated  with  the  pur- 
chase and  installation  of  an  absorption  chiller  is  also 
relatively  independent  of  solar  collector  area,  since 
the  selection  of  an  appropriate  absorption  chiller  is 
dictated  by  peak  building  cooling  load  for  all  but  the 
smallest  solar  collector  areas.  It  is  important  that 
these  costs  be  apportioned  appropriately  since  the 
method  compares  the  life-cycle  cost  difference 
between  conventional  and  solar  energy  systems  as  a 
function  of  the  area  of  the  solar  collector  array. 

When  determining  the  increased  capital  cost  of 
the  solar  system,  care  must  be  given  to  specifying  the 
type  of  auxiliary  system  to  be  used  and  to  taking  the 
appropriate  credits  for  any  reduction  in  capital  costs 
for  the  conventional  system.  Also,  the  costs  for  the 
auxiliary  energy  must  be  compatible  with  the 
auxiliary  system  being  used.  For  example,  most  solar 
systems  have  a complete  conventional  system  as  an 
auxiliary.  For  solar  heating  and  cooling  systems 
which  require  absorption  chilling,  the  auxiliary  can 
be  a boiler  for  the  absorption  chiller  or  an  auxiliary 
centrifugal  chiller.  In  the  first  case,  capital  cost  is 
low,  since  a credit  can  be  taken  for  the  centrifugal 
chiller  because  the  boiler  is  required  anyway  for  the 
heating  system.  The  fuel  costs,  however,  may  be 
high,  since  the  absorption  chiller  has  a COP  on  the 
order  of  .65.  In  the  second  case,  the  capital  costs  are 
higher,  since  credit  for  the  centrifugal  chiller  cannot 
be  taken.  However,  since  this  chiller  has  a COP  on 
the  order  of  4.2,  the  fuel  cost  may  be  considerably 
lower.  Note  that,  in  the  first  case,  heating  fuel  is  the 
auxiliary  fuel  while  in  the  second  case,  electrical 
power  is. 

Once  the  cost  data  have  been  determined,  the 
second  step  in  the  comparison  procedure  is  to  estab- 
lish a life-cycle  fuel  cost  for  conventional  energy 
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systems.  For  a conventional  system,  annual  fuel 
costs  can  be  estimated  as  follows: 

l.  Determine  the  annual  heating  and/or  cooling 
energy  requirement  for  conventional  equipment  as 
described  in  Chapter  2. 


Chapter  3.  Procedures  for  determining  life-cycle 
costs  for  various  solar  energy  systems  are  as  follows: 

1 . Determine  the  total  annual  auxiliary  energy  re- 
quirement for  a given  collector  area  and  system  type, 
Ola-  us'n8  the  method  described  in  Chapter  3. 


2.  Convert  the  cooling  energy  load  to  kilowatt 
hours  and  multiply  by  the  local  electrical  rate  to 
obtain  the  cost  of  cooling  energy  (note  that  the  COP 
of  the  conventional  electrical  chilling  device  must  be 
considered  in  making  this  conversion). 

3.  Divide  the  input  heating  energy  and/or  domes- 
tic hot  water  heating  energy  required  by  the  heating 
value  of  the  fuel  used  and  multiply  by  the  unit  fuel 
price  to  determine  the  cost  of  heating  (note  again 
that  in  determining  the  heat  energy  required,  the 
efficiency  of  the  boiler  or  furnace  must  be  con- 
sidered). 

4.  Add  results  of  steps  2 and  3 to  obtain  the  total 
annual  fuel  cost  for  the  conventional  system. 

5.  Convert  the  annual  fuel  cost  to  life-cycle  cost. 
This  value  provides  the  baseline  for  comparison  of 
solar  energy  systems. 

The  third  step  is  to  calculate  the  amount  of  aux- 
iliary fuel  or  electrical  energy  required  annually  by 
the  solar  system  being  considered  for  various  collec- 
tor array  sizes  based  on  the  methods  described  in 


2.  Determine  the  annual  energy  cost  for  auxiliary 
heat  similarly  to  conventional  systems.  For  heating 
and  cooling  systems,  assume  that  the  auxiliary 
energy  requirement  is  for  cooling  only.  Note  in  par- 
ticular that  if  centrifugal  or  other  electrical  chilling 
devices  are  to  be  used  to  augment  the  solar  cooling 
system,  the  appropriate  COPs  for  absorption  and 
electrical  refrigeration  equipment  must  be  used  in 
determining  the  auxiliary  energy  that  will  be  re- 
quired to  meet  the  cooling  demand  not  met  by  the 
solar  system.  If  a boiler  is  used  to  drive  the  absorp- 
tion chiller,  the  boiler  efficiency  must  be  used. 

3.  Convert  the  annual  auxiliary  fuel  cost  to  life- 
cycle  cost. 

The  total  life-cycle  cost  difference  between 
conventional  and  solar  heating  and  cooling  systems 
can  now  be  determined  for  various  collector  areas. 
For  a given  collector  area,  this  difference  is  simply 
the  life-cycle  capital  cost  difference  of  the  solar 
energy  system  plus  the  life-cycle  fuel  cost  for  the  aux- 
iliary system  minus  the  life-cycle  fuel  cost  of  the  con- 
ventional system.  Figure  4 is  the  hypothetical  plot  of 
this  cost  difference  for  varying  collector  areas. 
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Figure  4.  Net  life-cycle  cost  difference  vs.  collector  area. 
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Appendix  B provides  a detailed  example  which  fur- 
ther illustrates  application  of  this  technique. 

\ 

The  shape  of  the  curve  in  Figure  4 is  fairly  easy  to 
visualize.  For  very  small  collector  areas,  there  is  very 
little  fuel  savings  to  offset  the  fixed  costs  of  the 
system  such  as  the  controls,  added  piping,  and 
pumps,  etc.  As  collector  area  increases,  the  in- 
creased fuel  savings  start  to  offset  the  fixed  costs  of 
the  system,  and  the  cost  difference  between  solar  and 
conventional  costs  decreases.  As  more  collectors  are 
added,  a point  of  diminishing  returns  is  reached; 
added  collector  area  adds  very  little  additional  fuel 
savings.  Thus,  the  added  costs  of  collectors  and  tank 
volume  cannot  be  offset  by  the  fuel  savings  of  those 
collectors,  and  solar  system  cost  starts  to  increase 
with  respect  to  conventional  system  costs. 

When  the  cost  difference  for  the  solar  energy 
system  being  considered  has  been  plotted  in  the  form 
of  Figure  4.  the  optimum  collector  area  becomes 
obvious;  it  is  the  area  at  which  the  net  life-cycl?  cost 
difference  between  conventional  and  solar  heating 
and  cooling  systems  is  minimized.  Note  that  Figure  4 
concisely  summarizes  the  economics  of  solar 
systems.  The  points  on  the  curve  having  a positive 
cost  difference  indicate  that  solar  energy  is  not  cost 
effective.  Only  when  the  curve  dips  below  the  origin 
can  a solar  energy  system  be  economically  justified. 
Based  on  present  solar  energy  system  costs,  many 
applications  may  not  have  such  a negative  cost 
difference.  In  fact,  if  the  costs  of  collectors  and  tanks 
are  so  high  that  the  fuel  savings  cannot  offset  even 
those  costs  which  vary  with  collector  area,  the  curve 
in  Figure  4 will  increase  continuously. 

Based  on  the  methodology  presented  in  the  last 
three  chapters,  the  energy  and  life-cycle  cost  analysis 
of  the  solar  energy  system  can  be  made.  If,  based  on 
this  analysis,  the  solar  system  is  considered  to  be 
feasible,  then  a detailed  design  analysis  should  be 
made  using  the  computer-aided  method  described  in 
the  next  chapter. 


PERFORMING  COMPUTER-AIDED 
DESIGN  CALCULATIONS 

Introduction 

The  previous  chapter  described  a method  for  esti- 
mating graphically  the  performance  of  candidate 
solar  domestic  hot  water,  building  heating,  or  build- 
ing heating  and  cooling  systems.  This  method  is 


sufficiently  accurate  for  the  preliminary  energy  and 
life-cycle  cost  analysis  of  solar  energy  systems.  For 
most  projects,  however,  selecting  the  optimum  col- 
lector area  and  type  and  insuring  the  selection  of  the 
appropriate  collector  tilt  angle  and  storage  tank 
volume  may  require  a more  detailed  design  analysis. 
The  methods  described  in  this  chapter  are  based  on 
the  application  of  the  BLAST  program  developed  at 
CERL.  This  program  permits  hourly  simulation  of 
building  loads,  and  of  air-handling  systems  meeting 
the  various  loads  in  the  building,  and  simulation  of 
central  energy  plant  equipment  supplying  the  air- 
handlers,  including  solar  energy  systems.  Designers 
who  wish  to  use  this  program  should  obtain  the  users 
and  reference  manuals  which  describe  the  prepara- 
tion of  the  input  required  for  using  the  program  and 
the  algorithms  employed  in  performing  the  simula- 
tion. 

The  optimum  design  using  the  simulation  model  is 
again  determined  by  performing  the  economic 
analysis  as  described  in  Chapter  4.  The  simulation 
tools,  however,  allow  better  determination  of  system 
performance  and,  thus,  a more  accurate  measure  of 
q,  the  fraction  of  the  load  met  by  the  solar  system  for 
a given  collector  array.  Also,  other  variables  such  as 
collector  tilt  angle,  storage  tank  volume,  and  type  of 
collector  can  be  varied  to  maximize  q.  In  addition, 
the  BLAST  program  allows  for  the  analysis  to  be 
performed  using  more  accurate  values  of  the  build- 
ing energy  requirements,  Q^,  since  assumptions 
concerning  average  efficiency  and  COP  do  not  have 
to  be  made. 

Information  Required  to  Use 
the  BLAST  Program 

To  use  the  BLAST  program,  hourly  weather  data, 
including  solar  radiation  data  if  available,  must  be 
acquired  in  the  required  form  as  input  to  the  com- 
puter program.  Users  should  contact  CERL  to 
obtain  the  necessary  data  tapes  (CERL  will  arrange 
to  obtain  these  tapes  from  the  Air  Force’s  Air 
Weather  Service  or  from  the  NOAA).  These  tapes 
are  typically  referred  to  as  National  Weather  Service 
Series  280,  Solar  Data  Tapes  and  National  Weather 
Service  Series  1440,  Climatological  Data  Tapes. 
Obtaining  the  tapes  may  take  up  to  6 weeks;  there- 
fore, users  should  anticipate  use  of  the  program  and 
request  the  tapes  early  enough  so  they  will  be  availa- 
ble when  required. 

The  BLAST  program  requires  a description  of  the 
building  as  input.  The  users  manual  describes  the 
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exact  form  of  this  input.  The  program  uses  such 
information  as  the  type  of  wall  construction,  window- 
orientation,  schedule  for  occupancy,  and  electrical 
demands  in  determining  the  hourly  building  zone 
loads  and  hourly  total  load  of  the  building.  If  meas- 
ured data  are  available,  they  can  frequently  be  used 
to  check  the  building  loads  to  insure  that  parameters 
were  properly  defined  and  retlect  the  actual  build- 
ing. since  some  parameters,  such  as  infiltration  and 
occupancy,  are  hard  to  determine. 

The  air-handling  systems  used  in  the  building 
must  also  be  described.  In  this  case,  such  informa- 
tion as  the  zone  air  flow  rate,  type  of  coil,  type  of  air- 
handling system,  and  type  of  control  system  used  are 
specified  by  the  user.  The  output  from  the  air  distri- 
bution system  simulations  is  the  hourly  energy- 
demand  on  the  boilers  and  chillers  serving  the  air 
handlers  in  the  building.  These  hourly  demands 
form  the  input  to  the  central  plant  simulation  pro- 
gram which  determines  the  expected  performance  of 
conventional  and  solar  energy  system  components. 
In  performing  the  central  energy  plant  simulation, 
the  user  can  specify  the  performance  of  the  solar 
collector,  the  collector  area,  and  the  storage  tank 
volume  being  used  in  the  system.  The  type  of  per- 
formance data  required  for  solar  collectors  is  identi- 
cal to  that  described  in  National  Bureau  of  Stand- 
ards publication  NBSIR-74-635.10  The  slope  and 
intercept  of  solar  collector  efficiency  curves  are  the 
input  parameters  required. 

Note  that  in  the  analysis  of  the  solar  energy 
systems,  the  BLAST  program  should  first  be  used  to 
minimize  the  energy  required  by  each  zone  by  study- 
ing such  effects  as  added  insulation  and  changes  in 
window  area  and  orientation.  This  should  be 
followed  by  a comparison  of  various  air-handling 
systems  and  control  schemes  to  minimize  the  energy 
demanded  by  the  boilers  and  chillers,  assuming 
these  design  changes  can  be  made  to  the  building. 
Once  this  analysis  has  been  performed,  various 
collector  systems  must  be  simulated  to  optimize  the 
performance  of  the  candidate  systems  and  to  estab- 
lish the  optimum  system  configuration. 

Once  the  load  profiles  for  the  boilers  and  chillers 
have  been  determined,  the  optimization  of  the  solar 
energy  system  can  be  initiated. 


,0J.  S.  Hill  and  T.  Kusuda.  Methods  of  Testing  for  Rating  Solar 
Collectors  Based  on  Thermal  Performance.  NBSIR-74-635 
(National  Bureau  of  Standards.  1474). 


Setting  Up  the  Solar  Energy 
Detailed  Design  Study 

The  computer-aided  design  should  not  be  initiated 
until  the  graphical  life-cycle  cost  analysis  described 
in  the  previous  chapter  has  been  completed,  since 
the  results  of  that  study  provide  the  starting  point  for 
the  computer-aided  study.  If  the  life-cycle  cost  study 
has  been  completed,  but  the  loads  determined  from 
the  BLAST  program  are  considerably  different  than 
those  used  in  the  life-cycle  cost  analysis,  the  analysis 
should  be  repeated.  The  graphical  method  serves  to 
bracket  the  area  of  consideration. 

In  order  to  optimize  the  design,  certain  variables 
not  previously  considered  should  be  varied  to  estab- 
lish their  optimum  value,  as  discussed  in  the  follow- 
ing section.  The  starting  point  should  be  the  solar 
system  obtained  from  the  universal  curve  analysis 
which  set  the  optimum  collector  area.  The  collector 
array  should  be  set  with  an  azimuth  angle  of  180 
degrees  (i.e.,  facing  due  south)  and  tilted  at  the  opti- 
mum angle  as  described  in  Chapter  2.  The  tank 
volume  should  be  set  at  16  lb/sq  ft  of  collector 
(78  kg/m2). 

Methodology  for  Optimizing 
the  Solar  System  Design 

The  optimization  proceeds  in  the  following  steps: 

1.  First,  several  simulations  are  run  using  the 
BLAST  solar  simulation  package  with  different  tilt 
and  azimuth  angles  for  the  collector  array.  Figure  5 
illustrates  how  the  output  from  such  a parametric 
study  can  be  plotted  to  establish  the  optimum  tilt  for 
a given  azimuth.  A similar  plot  of  system  perform- 
ance versus  azimuth  angle  will  indicate  the  optimum 
azimuth. 

2.  The  next  step  is  to  vary  the  collector  area  while 
keeping  the  collector-area-to-tank-volume  ratio  con- 
stant (start  with  16  lb/sq  ft  of  collector  [78  kg/m2]) 
and  repeat  the  economic  analysis  described  in  Chap- 
ter 4,  obtaining  a curve  similar  to  Figure  4.  The 
lowest  dip  in  the  curve  is  the  optimum  collector  area 
for  this  collector-area-to-tank-volume  ratio.  (Note 
that  this  is  a repeat  of  the  universal  curve  analysis 
using  the  more  precise  computer  analysis  tools.) 

3.  The  third  step  is  to  repeat  the  previous  step  for 
different  collector-area-to-storage-tank-volume  ra- 
tios. All  life-cycle  cost  curves  can  be  plotted  on  the 
same  graph.  Sufficient  ratios  should  be  examined  to 
determine  the  curve  with  the  lowest  dip.  If  no  other 
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Figure  5.  System  performance  vs.  collector  slope. 


information  is  available,  8 and  32  Ib/sq  ft  (39  and 
156  kg/m2)  are  good  second  and  third  tries.  Subse- 
quent tries  are  based  on  the  results  of  these  simula- 
tions. The  curve  with  the  lowest  dip  gives  the  opti- 
mum tank-volume-to-collector-area  ratio,  while  the 
point  at  which  the  dip  occurs  establishes  the  opti- 
mum collector  area.  Although  the  better  the  tank 
insulation,  the  larger  the  optimum  tank  volume 
will  be.  high  tank  insulation  (especially  for  buried 
tanks)  is  unfortunately  difficult  to  achieve  in  prac- 
tice. Therefore,  tank-volume-to-collector-area  ratios 
which  greatly  exceed  16  Ib/sq  ft  (78  kg/m2)  should 
be  examined  critically  to  insure  the  specified  tank 
insulation  values  can  be  obtained. 

4.  The  final  step  is  to  vary  the  type  of  collector  to 
see  what  effects  this  has  on  system  life-cycle  costs. 
The  optimum  tank-volume-to-collector-area  ratio 
found  in  Step  3 should  not  be  significantly  affected 
by  the  type  of  collector.  The  life-cycle  cost  difference 
versus  collector  area  curves  should  be  constructed  as 
before  by  simulating  several  different  areas  for  each 
collector  type  being  considered. 

Note  that  the  only  differences  between  the  graphi- 
cal and  computer-aided  methods  for  determining 
system  performance  and  optimal  collector  configura- 


tion are  that  the  simulation  model  provides  a more 
accurate  estimate  of  the  performance  of  a candidate 
solar  energy  system  and  permits  a more  refined  study 
of  the  secondary  effects  of  collector  tilt,  azimuth 
angle,  and  storage  tank  capacity. 

6 PRACTICAL  CONSIDERATIONS 

Discussions  in  the  previous  chapters  focused  on 
estimating  performance  of  solar  collector  systems. 
This  chapter  deals  with  some  practical  considera- 
tions involved  in  the  final  design,  installation,  and 
operation  of  solar  energy  systems  and  is  based  on 
experience  gained  in  the  construction  and  operation 
of  the  CERL  solar  facility. 

Controls 

The  controls  of  solar  energy  systems  typically  do 
not  require  new  control  technology.  However,  special 
care  must  be  exercised  in  the  control  of  the  collector 
array.  The  most  convenient  way  to  control  the  collec- 
tor pump  which  determines  when  the  collector  is  or 
is  not  collecting  energy  is  by  sensing  the  temperature 
difference  between  the  collector  plate  surface  and 
the  storage  tank.  When  a positive  difference  exists 
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(i.e.,  the  collector  is  hotter  than  the  storage  tank), 
useful  energy  can  be  collected  and  the  collector 
system  should  be  energized.  However,  any  time  the 
collector  plate  is  cooler  than  the  storage  tank,  the 
system  will  dissipate  or  waste  heat  through  the  col- 
lector array  if  it  is  energized.  Thus,  care  must  be 
taken  to  insure  that  when  a negative  difference  exists 
between  the  temperature  of  the  collector  absorber 
plate  and  the  temperature  in  the  storage  tank,  the 
system  is  turned  off.  At  the  same  time,  minimizing 
the  difference  between  those  two  temperatures  when 
the  system  is  to  be  turned  on  is  desirable.  For 
example,  it  might  be  appropriate  to  establish  a 10°F 
(b°C)  difference  between  the  collector  plate  and  the 
storage  tank  for  energizing  the  system  and  a 5°F 
(3°C)  difference  for  de-energizing  it. 

It  is  also  important  that  the  control  system  main- 
tain its  linearity  over  the  entire  region  of  control.  For 
example,  turning  the  collector  loop  on  may  be  re- 
quired when  the  collector  is  185°F  (85°C)  and  the 
tank  is  I80°F  (82°C)  during  the  summer.  During  the 
winter,  the  temperature  may  be  I25°F  (52°C)  on  the 
collector  surface  and  120°F  (49°C)  in  the  tank;  thus, 
the  measurement  of  temperature  at  the  plate  and  in 
the  tank  and  the  subsequent  sensing  of  the  tempera- 
ture difference  between  the  two  require  linear  con- 
trol over  a broad  temperature  range. 

Piping  of  the  Auxiliary  Energy  Supply 

For  practical  solar  energy  systems,  some  auxiliary 
energy  supply  will  be  required.  Since  there  are  many 
possible  configurations,  individual  cases  must  be 
examined  to  determine  where  auxiliary  energy  will 
be  derived.  However,  any  auxiliary  supply  boiler  or 
other  device  should  be  piped  parallel  with  the  solar 
energy  solar  tank.  Figure  6 illustrates  the  correct  and 
incorrect  methods  for  piping  an  auxiliary  energy 
supply.  If  the  auxiliary  supply  is  piped  in  series  with 
the  tank,  under  conditions  of  low  tank  temperature 
the  auxiliary  supply  will  not  only  be  required  to  meet 
the  heat  energy  demands  of  the  building  but  will  also 
begin  to  add  energy  to  the  tank.  This  is  particularly 
undesirable  because  it  will  require  an  unnecessarily 
large  auxiliary  heater  capable  of  meeting  the  instan- 
taneous demand  while  at  the  same  time  supplying  an 
added  demand  by  delivering  energy  to  the  storage 
tank.  It  is  also  undesirable  because  the  stored  energy 
should  be  derived  only  from  the  solar  energy  system; 
application  of  auxiliary  heat  to  the  stored  energy 
raises  the  tank  temperature,  reducing  collector  effi- 
ciency and  increasing  tank  heat  losses. 


Use  of  Two-Position  Valves 

Another  consideration  in  the  practical  application 
of  solar  energy  is  that  free  convection  be  prevented 
from  dissipating  solar  heat  through  the  collector 
bank  when  the  solar  energy  system  is  not  in  opera- 
tion. The  application  of  two-position  valves  which 
close  when  the  collector  system  is  de-energized  is 
recommended  to  prevent  this. 

Venting  the  Collector  Loop 

Particular  attention  must  be  given  to  the  possi- 
bility of  air  accumulating  in  the  collector  array.  A 
frequent  problem  encountered  in  the  past  has  been 
accumulation  of  pockets  of  air  at  the  top  of  the  col- 
lector bank.  These  air  pockets  prevent  the  fluid  from 
passing  over  the  collector  absorber  plate,  thus  reduc- 
ing the  area  over  which  solar  energy  is  collected. 
Such  "vapor  locking"  can  be  prevented  by  generous 
use  of  automatic  or  manual  vents  in  the  collector 
system  and  by  insuring  that  sharp  bends  or  other 
piping  arrangements  which  can  lead  to  air  being 
trapped  are  avoided.  Air  separators  in  the  collector 
loop  are  also  recommended. 

Piping  of  Hot  Fluids 

Additional  caution  should  be  exercised  in  posi- 
tioning pumps  in  the  system,  particularly  where  the 
storage  tank  is  to  be  maintained  at  or  near  atmo- 
spheric pressure.  A potential  cause  of  system  failure 
is  the  cavitation  of  pumps  used  to  pump  hot  fluids. 
Piping  considerations  should  include  assurances  that 
both  the  collector  loop  and  storage  tank  are  posi- 
tioned such  that  net  positive  suction  exists  at  the 
pump  at  all  times  to  prevent  flashing  when  the 
pumps  are  energized.  Self-priming  pumps  should  be 
used  whenever  appropriate. 

Storage  Tank  Insulation 

Practical  experience  indicates  that  the  heat  loss 
from  insulated  tanks  is  frequently  higher  than  pre- 
dicted. For  buried  tanks  in  particular,  special  care 
must  be  taken  to  insure  that  an  impenetrable  vapor 
barrier  surrounds  the  insulation  on  the  tank.  Con- 
servative buried  tank  volumes  should  be  considered 
to  minimize  the  heat  loss.  This  consideration  is  less 
important  for  tanks  located  inside  the  building. 
However,  such  tanks  occupy  valuable  floor  space 
whose  worth  should  be  considered  in  analyzing  the 
cost  of  the  tank.  Also,  in  applications  where  cooling 
is  considered,  the  location  of  tanks  inside  the  build- 
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Figure  6.  Piping  the  auxiliary  energy  supply. 


ing  must  be  considered  so  that  they  do  not  contribute 
to  the  overall  cooling  load  in  the  summer. 

7 CONCLUSIONS 

The  following  conclusions  can  be  drawn  concern- 
ing the  methods  presented  in  this  report: 

1.  The  graphical  method  of  evaluating  solar 
energy  systems  is  based  on  established  technology 
and  calculation  methods,  and  provides  sufficient 
accuracy  for  use  in  feasibility  studies. 

2.  The  life-cycle  cost  analysis  method  includes  all 
cost  factors  which  must  be  considered  in  evaluating 
solar  energy  systems,  and  interfaces  directly  with  the 
OCE  life-cycle  cost  analysis  instructions. 

3.  The  computer  simulation  method  is  based  on 
established  calculation  methods  and  can  be  used  in 
preparing  design  instructions  and  evaluating  final 
designs. 

4.  The  input  data  required  for  both  the  graphical 


and  computer  simulation  methods  can  be  readily 
obtained  by  District  personnel. 

8 FUTURE  PLANS 

Both  the  graphical  and  computer  procedures  for 
making  a solar  system  energy  and  life-cycle  cost 
analysis  will  be  field  tested  in  fiscal  year  1977  (FY  77) 
at  a Corps  of  Engineers  District  office  which  is  cur- 
rently performing  a solar  system  study.  The  field  test 
will  determine  if  the  procedures  are  clear  and  user- 
oriented  and  whether  the  input  data  are  available. 
Based  on  the  results  of  the  field  test,  the  procedures 
will  be  revised  to  make  them  more  understandable 
and  easier  to  use. 

The  graphical  and  computer  simulation  perform- 
ance evaluation  tools  will  be  validated  using  actual 
data  from  solar  demonstrations  at  CERL  and  Army 
installations  during  FY  77  and  FY  78. 

During  FT  78  the  revised  procedures  will  be 
implemented  by  Corps  Districts,  and  training  ses- 
sions in  their  use  will  be  given. 
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APPENDIX  A: 


SOLAR  RADIATION  DATA* 


•From  Climatic  Atlas  of  the  United  States  (National  Oceanic 
and  Atmospheric  Administration.  I%8). 


Figure  A2.  Mean  daily  solar  radiation  tor  January. 
Conversion  factor:  1 langley  = 3.69  Btu/sq  ft  or  41 .84  kj 


I 


radiation  for  March 


Figure  A8.  Mean  daily  solar  radiation 
Conversion  factor:  1 laneley  = 3.69  Btu/sq  ft 
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Figure  A9.  Mean  daily  solar  radiation  for  August. 
Conversion  factor:  1 langley  = 3.69Btu/sq  ft  or  41 .84  kJ/m2. 


Figure  A10.  Mean  daily  solar  radiation  for  September. 
Conversion  factor:  1 langley  = 3.69  Btu/sq  ft  or  41.84  kJ/m 


Figure  All.  Mean  daily  solar  radiation  for  October. 
Conversion  factor:  1 langley  = 3.69  Btu/sq  ftor41.84kJ/mJ. 


Figure  A12.  Mean  daily  solar  radiation  for  November. 
Conversion  factor:  I langlev  = 3.69  Btu/sq  ft  or  41 .84  kJ/m!. 


APPENDIX  B: 

EXAMPLE  OF  SOLAR  SYSTEM 
LIFE-CYCLE  COST  ANALYSIS 

A design  engineer  wants  to  determine  the  collector 
area  which  is  economically  the  most  feasible  for  solar 
heating  and  cooling  of  a new  administration  building 
similar  to  the  one  used  for  Examples  2 and  3 in 
Chapter  3.  For  the  purpose  of  this  example,  he/she 
is  to  compare  the  costs  of  a conventional  building 
with  oil  heating  and  electric  centrifugal  cooling  to  a 
building  heated  and  cooled  by  a solar  energy  system. 
In  the  solar  case,  auxiliary  heating  is  supplied  by  an 
oil-fired  boiler  and  auxiliary  cooling  is  provided  by 
an  electric  centrifugal  chiller. 

Input  Data: 

From  Examples  2 and  3 the  annual  loads  are 
Building  Heating  Load  108  x 10'  Btu/yr* 
Building  Cooling  Load  112x10'  Btu/yr 
The  cost  of  energy  is  assumed  to  be 
Price  of  Electricity  $. 027/kWh 

Price  of  Oil  $.5/gal**t 

Heating  and  Cooling 
Step  I 

Following  the  procedure  on  pp  18-20  of  Chapter  4, 
the  capital  costs  of  the  solar  system  components 
must  be  determined.  Since  in  this  case  the  absorp- 
tion chiller  is  redundant,  its  entire  price  is  con- 
sidered. Table  B1  summarizes  estimated  costs. 

Since  the  building  of  Examples  2 and  3 contains  a 
385-ton  cooling  unit,  the  cost  of  an  absorption  chiller 
($470/ton)  is  $1.81  x 10s.  This  is  the  value  listed  in 
the  table. 


•SI  conversion  factor:  1 Btu  - 1 .055  kl.  (To  avoid  confusion,  SI 
equivalents  are  not  given  in  this  example;  instead,  an  SI  conver- 
sion factor  is  given  the  first  time  a particular  unit  is  used.) 

••SI  conversion  factor:  1 gal  = 3.785  ,. 
t'This  price  is  used  for  example  purposes  only;  it  may  be  some- 
what higher  than  the  actual  fuel  cost. 


Table  Bl 

Estimated  Costs  for  Heating  and  Cooling  Example 


Fixed  Costa,  $ 

Coat/ iq  ft* 
of  Collector,  $ 

Absorption  Cooling 

System 

1.81  x I04 

Controls 
Heat  Exchanger 

3 » I0> 

.5 

fluid 

.18 

Storage  Tank 

1.0 

Plumbing 

.75 

Pumps 

.75 

Collector 

5.0 

Labor 

1.0 

SI. 84*  l()5 

S9.I8AC 

•SI  conversion  factor:  I sq  ft  = 0.0029  m1. 

The  cost  of  control  sensors  and  valves  is  nearly 
independent  of  system  size  and  is  estimated  at 
roughly  $3000. 

The  remainder  of  the  items  in  the  table  are 
dependent  on  collector  area  <AC).  These  values  are 
meant  as  estimates  for  the  purposes  of  the  example 
only. 

Using  the  OCE  life-cycle  costing  instructions,  the 
life-cycle  initial  costs  can  be  expressed  as 

CC  = Y(1.84xi05  + 9.18Ac) 

where  Y = the  initial  cost  multiplying  factor.11  For 
this  example,  Y is  taken  to  be  2.05.  assuming  an 
interest  rate  of  6. 5 percent  and  a 25-year  system  life. 

Step  2 

Once  the  cost  data  have  been  determined,  the 
second  step  in  the  comparison  procedure  is  the 
establishment  of  a life-cycle  fuel  cost  for  conven- 
tional energy  systems.  As  described  in  Chapter  4, 
conventional  system  annual  fuel  costs  are  estimated 
as  follows: 

1.  Determine  the  annual  heating  and/or  cooling 
energy  requirement  for  conventional  equipment  as 
described  in  Chapter  2. 

For  this  example  the  heating  load  is  108  x 10*  Btu/ 
yr  and  the  cooling  load  is  112  x 10'  Btu/yr,  as  was 
shown  in  Examples  2 and  3 of  Chapter  3. 

11 OCE  Life  Cycle  Costing  Instructions  (Department  of  the 
Army.  May  1971). 
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2.  Convert  the  cooling  energy  load  to  kilowatt  Step  .? 
hours  and  multiply  by  the  local  electrical  rate  to 

obtain  the  cost  of  cooling  energy  (note  that  the  The  third  step  is  to  calculate  the  amount  of  aux- 

coefficient  of  performance  of  the  conventional  elec-  iliary  fuel  or  electrical  energy  required  annually  by 

trical  chilling  device  must  be  considered  in  making  the  solar  system  being  considered  for  various  collec- 

this  conversion).  tor  array  sizes  based  on  the  methods  described  in 

Chapter  3.  Procedures  for  determining  life-cycle 
For  a cooling  load  of  1 12  x 10*  Btu/yr  and  a chiller  costs  for  various  solar  energy  systems  are  as  follows: 
COP  of  4.  the  cost  of  cooling  energy  (at  $. 027/kWh) 

is:  1.  Determine  the  total  annual  auxiliary  energy 

requirement  for  a given  system,  Ola-  by  the  method 
<11?  x tn*  Rm/vr)  (1/4)  f kWh  \fS  027\  prescribed  in  Chapter  3.  Following  Example  2 of 

- ' \3412  Btu/  \ kWh  / Chapter  3,  Ql  is  280  x 10*  Btu/yr  and  Qc  is  5.6  * 105 

Btu/sq  ft/yr.  For  a collector  area,  Ac,  of  10s  sq  ft  (a 
= S2.2  x 104/yr  guess),  Eq  1 may  be  used  to  find  Ps 

3.  Divide  the  input  heating  energy  and/or  domes- 
tic hot  water  heating  required  by  the  heating  value  of 
the  fuel  used  and  multiply  by  the  unit  fuel  price  to 
determine  the  cost  of  heating  (note  again  that  in 
determining  the  heat  energy  required,  the  efficiency 
of  the  boiler  or  furnace  must  be  considered). 


Ps  = 2.0 

For  this  Ps,  the  universal  curve  for  heating  and  cool- 
ing (Figure  2)  gives  a q of  .49.  Hence. 

Ola  = U-e)0L  = d--49)0. 


For  a boiler  which  is  75  percent  efficient,  using  a Of  a = 143  x 10s  Btu/yr 

heating  value  of  150,000  Btu/gal  for  oil,  the  cost  of 

heating  (at  S. 5/gal  for  oil)  is  2.  Determine  the  annual  energy  cost,  Fs.  for  aux- 

iliary fuel.  (All  auxiliary  energy  is  assumed  to  be 
(108  Btu/  r)  ( - V*'5  V * ga,_  \ used  for  cooling.)  Since  a centrifugal  chiller  (COP  of 

y ' \.75/\gal  /\1.5x  10s  Btu/  4)  supplies  the  cooling  backups,  the  auxiliary  fuel 

requirement  may  be  determined  as  follows: 

= $4.8  x lOVyr 

Fs  = $1 .84  x 104/yr 

4.  Add  results  of  steps  2 and  3 to  obtain  the  total 

annual  fuel  cost,  Fc,  for  the  conventional  system:  The  factor  of  .65  comes  from  the  fact  that  Ql 

(280  x 1 08)  was  computed  on  the  basis  of  absorption 
Fc  = 2.2  x 104  + 4.8  x 104  = $7.0  x 104/yr  chiller  cooling.  Now  that  the  backup  system  is 

assumed  to  be  conventional,  this  factor  is  regained. 

5.  Convert  the  annual  fuel  cost  to  life-cycle  oper- 
ating cost.  This  value  provides  the  baseline  for  com-  3.  Convert  the  annual  fuel  cost  to  life-cycle  oper- 

parison  of  solar  energy  systems.  ating  cost  LOs.  Once  again 


The  life-cycle  operating  cost  for  fuel  for  the  con- 
ventional system  (LOe)  is  given  by 

LOc  = M Fc 


LOs  = MFS  = $0.78  x 106 

Here  the  value  of  M is  the  same  as  for  the  conven- 
tional system. 


where  M = a multiplier  which  gives  life-cycle  fuel  Step  4 
costs.  For  oil  and  a 25-year  facility  life,  M = 42.4. 

Hence,  The  total  life-cycle  cost  difference  (LCC)  between 

it  < nventional  and  solar  heating  and  cooling 
LOc  = (42.4)  ($7.0  x |04)  systtc:  .s  given  by 


LOc  = $2.97  x 10* 


LCC=CC  + LOs-LOc 
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The  solar  energy  system  represents  a savings  over 
the  conventional  system  only  if  LCC  is  negative. 

From  Step  1 , for  an  Ac  of  1 x 105  sq  ft 

CC  = 2.0S|$1.84x  105  + $9.18(1  x 105)| 

CC  = $2.26  x 10‘ 

Using  LOs  and  LOc  from  Steps  2 and  3,  LCC  can 
now  be  computed. 

LCC  = $2.76  x 10‘  + .78  x 10‘  -2.97  x 10* 

LCC  = .07  x 10‘ 

for  a collector  area  of  100.000  sq  ft. 

A life-cycle  cost  difference  of  this  magnitude  is 
negligible  when  compared  to  the  capital  or  fuel  life- 
cycle  costs.  Hence,  to  within  the  accuracy  of  the 
method,  a solar  energy  system  of  100,000  sq  ft  is  (for 
this  case)  equivalent  to  a conventional  system.  Other 
collector  areas  may  yield  greater  (or  smaller)  differ- 
ences. 

If  the  previous  calculation  is  repeated  for  an  Ac  of 
1.5  x 105  sq  ft,  it  is  found  that 

Ps=3 

The  universal  curve  gives  a (/  of  .67  for  this  Ps;  thus, 
Ola  = 92  x 10*  Btu/yr 
In  the  same  way 

Fs  =$1.18  x I04 

This  implies 

LOs  = $.5  x 10‘ 
for  an  Ac  of  1.5  x 10s  sq  ft 

CC  =$2.05|1.84  x 10s  + 9. 18(1. 5 x 105))  =$3.2x  10* 
Using  the  same  L Oc  as  before 

LCC  = 3.2  x 10‘ + .50  x IO»-2.97x  10‘ 

LCC  = $.73  x I0‘ 

for  a collector  area  of  150,000  sq  ft.  The  larger  col- 
lector area  has  resulted  in  a higher  life-cycle  cost. 


If  a collector  area  of  75,000  sq  ft  is  considered, 
LCC  is  computed  to  be  — $2.3  x 10*.  Since  a cost- 
effective  area  has  been  found,  a more  exact  analysis 
of  the  problem  is  indicated.  Hence,  the  computer 
simulation  program  should  probably  be  consulted. 

Heating  Only 

Experience  indicates  that  a high  percentage  of  the 
auxiliary  fuel  requirement  is  normally  used  for  cool- 
ing. This  seems  to  suggest  that  a heating-only  system 
might  (from  an  economic  standpoint)  be  more  feasi- 
ble. An  example  of  a heating-only  application 
follows. 

Step  I 

Determine  capital  costs.  Table  B2  summarizes 
estimated  costs.  For  this  case,  the  absorption  chiller 
cost  is  saved,  and  the  life-cycle  capital  (initial)  cost, 
CC,  is  given  by 

CC  = Y(3x  10J  + 9.81  Ac) 


Table  B2 

Estimated  Cotta  for  Heating-Only  Example 


Fixed  Co«t,  $ 

Coat/aq  ft 
of  Collector,  $ 

Controls 

3 x 101 

Heat  Exchanger 

.5 

Fluid 

.18 

Tank 

1.0 

Plumbing 

.75 

Pumps 

.75 

Collector 

5.0 

Labor 

1.0 

3x  KF 

S9.18AC 

Step  2 


The  computation  of  conventional  fuel  costs  is 
based  on  the  heating  load  alone.  For  a building  heat- 
ing load  of  108  x 101  Btu/yr,  with  a .75  boiler  effi- 
ciency, 

Fc  = (108  x 10*  Btu/yr)  (1/.75)  ($.5/gal) 

(1  gal/1. 5x  10s  Btu) 

Fc  = $4.8  x 104/yr 

Hence 


LO„  = $2,035  x 10‘ 
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Table  B3 

Summary  of  Procedure  Used  to  Determine 


Month 

Q^,  Btu/mo 

Qc,  Btu/sq  ft/mo 

P. 

Q 

QL£,  Btu/mo 

Oct 

2.6  x 10* 

4.7  x |0< 

9.04 

.99 

2.6  x 10* 

Nov 

8.4  x 10* 

3.8  x 10* 

2.26 

.59 

S.Ox  10* 

Dec 

I2.9x  10* 

3.5  x 10* 

1.36 

.39 

S.Ox  10* 

Jan 

20.5  x 10* 

4.0  x 10* 

.98 

.29 

5.9  x 10* 

Feb 

19.1  x 10" 

4.5  x 10* 

1.18 

.35 

6.7x  10* 

Mar 

I6.9x  10* 

4.9  x 10* 

1.45 

41 

6.9x  10* 

Apr 

lJ.Sx  10* 

4.4  x 10* 

1.63 

.45 

6.1  x |0« 

May 

9.4  x 10* 

4.9x  10* 

2.67 

.65 

6.1x10* 

Jun 

4.6  x 10* 

4.8  x 10* 

5.22 

.97 

4.5x  10* 

48.8  x 10*  Btu 

Step  3 

To  compute  the  cost  of  auxiliary  fuel  for  the  solar 
system,  Qla  must  be  determined  from  monthly  Ql 
and  Qc  data.  The  values  for  Ql  and  Qc  for  the  heat- 
ing months  are  given  in  Example  3,  Chapter  3.  This 
information  (along  with  an  assumed  collector  area) 
allows  calculation  of  a Ps  for  each  month.  The  uni- 
versal curve  for  heating  (Figure  3)  then  gives  the  q 
for  each  Ps.  Finally  a monthly  QlS  can  be  figured 
(pQ i)  and  summed  to  get  a yearly  Ols- 

Table  B3  summarizes  the  procedure  using  the  Ql 
and  Qc  from  Example  3 of  Chapter  3,  with  a collec- 
tor area  of  5 x 104.  The  Qla  can  be  determined  from 

QLa  = 0L  - 0LS  = 108  x 10*  - 49  x 10* 

Qla  = 59  x 10*  Btu/yr 

The  solar  auxiliary  fuel  cost,  Fs,  is  given  (assuming  a 
boiler  efficiency  of  .75)  by 

Fs  = (59  x 108  Btu/yr)  (l/.75)($.5/gal) 

(1  gal/1.5  x 10*  Btu) 


LCC=CC  + LOs-LOc 
for  an  Ac  of  5 x ft)4  Sq  ft 
CC  = 2.05  (3xl0J  + 9.l8Ac) 

= 2.05|3  x 103  + 9. 18(5  x 104))  = S9.47  x 10s 

Thus, 

LCC  = S9.47  x 10‘  + $1. 13  x 10‘  — 2.04  x |0‘ 
LCC  = $.037  x 106 

for  a collector  area  of  50,000  sq  ft.  Once  again,  the 
value  is  small  when  compared  to  the  capital  and  fuel 
life-cycle  costs.  The  price  of  the  solar  energy  system 
compares  favorably  to  the  cost  of  the  conventional 
system . 

Table  B4  shows  results  of  the  analysis  for  other 
collector  areas.  The  solar  system  of  maximum  cost 
benefit  is  seen  to  occur  at  a collector  area  of  roughly 
25,000  sq  ft.  At  this  point  the  computer  simulation 
program  should  be  consulted  for  a more  exact 
analysis. 


Fs  = $2.7  x 104/yr 


Table  B4 


Hence, 


LOs  = MFX=  (42.4)  (2.7  xlO4) 
LOs  = $1 .13  x 10* 


Step  4 

Determine  LCC,  where 


Life-Cycle  Costs  for  Various  Collector  Areas 


Collector  Area  Isq  ft) 

LCC  (S) 

0 

+ .6x  10* 

10.000 

— 3.6  x |0* 

25.000 

-5.3  x 10* 

40.000 

-2.1  x u>* 

50.000 

+3.7x  HI* 

75.000 

+3+I.9X  |0» 
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